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We derive a quantization formula of Bohr-Sommerfeld type for computing quasinor- 
mal frequencies for scalar perturbations in an AdS black hole in the limit of large scalar 
mass or spatial momentum. We then apply the formula to find poles in retarded Green 
functions of boundary CFTs on IR^''^"^ and IR x S'^~^ . We find that when the boundary 
theory is perturbed by an operator of dimension A ^ 1, the relaxation time back to equi- 
librium is given at zero momentum by ^ ^ . Turning on a large spatial momentum 
can significantly increase it. For a generic scalar operator in a CFT on 1R^''^"\ there exists 
a sequence of poles near the lightcone whose imaginary part scales with momentum as 

d-2 

p d+2 in the large momentum limit. For a CFT on a sphere S we show that the theory 
possesses a large number of long-lived quasiparticles whose imaginary part is exponentially 
small in momentum. 
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1. Introduction 

Black hole quasinormal modes are damped oscillations which describe the evolution 
of small classical perturbations in a black hole background. These perturbations die off by 
falling into the horizon of the black hole or, in the case of an asymptotically flat spacetime, 
by escaping to infinity. The frequency of the oscillations therefore is not real but has 
an imaginary part describing their damping rate. Study of such small perturbations are 
important for understanding the stability of a black hole and have important applications 
for gravitational wave astronomy (see e.g. for a review). 

In the AdS/CFT correspondence [0,^3 5 a finite temperature boundary gauge 
theory is described by a black hole in AdS 0, study of quasinormal modes for an AdS 
black hole is important for understanding dynamics of strongly coupled gauge theories 
at finite temperature More explicitly, black hole quasinormal modes correspond 

in the dual gauge theory to poles of thermal retarded Green functions in the complex 
frequency plane Thus the decay of the black hole quasinormal modes describes in 

the gauge theory how small perturbations of the thermal state relax back to equilibrium. In 
particular, when the imaginary part of a quasi-normal frequency is smaller than the inverse 
temperature and the spacing between the neighboring modes, it describes a quasiparticle 
in the boundary gauge theory. 

Singularity 

Boundary Hor i zon Boundary 

Singularity 

Fig. 1: Pole structure in the complex frequency {u) plane of the thermal Wight- 
man function of a typical operator in the boundary theory. The poles in the lower 
half plane are the quasi-normal poles for the retarded Green function. Those in 
the upper half plane are reflections of the ones in the lower half plane with re- 
spect to the real axis as follows from the standard relation between a retarded and 
Wightman function. Each pole line extends to infinity. The analytic behavior of 
the Wightman function in different regions of the frequency plane reflects the bulk 
geometry in different regions of the black hole spacetime. For example the analytic 
behavior near a; = reflects the geometry near the black hole horizon. 
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It has been argued in that the analytic structure of thermal correlation func- 

tions of the boundary theory encodes information about the region of spacetime behind 
the horizon and in particular the black hole singularities. Understanding the structure of 
quasinormal polesll was an important element in the analysis. For example, the thermal 
Wightman two-point function of a typical operator in the boundary theory has the pole 
structure shown in fig. |l|. The four lines of poles divide the asymptotic frequency plane 
into four regions. It was pointed out in [|11| that the analytic behavior of the Wightman 
function in different regions of the frequency plane reflects the bulk geometry in different 
regions of the black hole spacetime, i.e. there exists a mapping between the complex- 
frequency plane of the boundary gauge theory and the bulk geometry. Such a map can be 
considered as a generalization to regions inside the horizon of the standard IR/UV connec- 
tion. The regions near the imaginary and real infinities in the frequency plane are mapped 
to the regions near the black hole singularities and the boundaries respectively, as labeled 
in fig. In particular, the curvature divergence at the singularities is reflected in the 
divergences of derivatives of the YM response functions at large imaginary frequencie^. 

Motivated from understanding the pole structure in fig. |^ at general spatial momenta, 
in this paper we derivei a general quantization formula of Bohr-Sommerfeld type for de- 
termining the quasinormal frequencies of scalar perturbations of an AdS black hole in the 
limit of large scalar mass or large spatial momentum. The basic idea for our approxima- 
tion is as follows. We first rewrite the Laplace equation for the scalar field in terms of a 
Schrodinger equation using the tortoise coordinate of the black hole geometry. We then 
simplify the Schrodinger equation using a WKB approximation which can be considered as 
a generalization (to finite or large masses) of the geometric optics approximation for mass- 
less modes. Quasinormal modes lie on anti-Stokes lines in the complex frequency plane 
where a certain sub dominant correction becomes comparable to the leading WKB result 
and can be determined by requiring the two contributions to cancel each other exactly. 



^ Given the one-to-one correspondence between quasi-normal modes and poles of retarded 
correlators in the boundary, we will refer to poles in gauge theory correlators as quasi-normal 
poles to emphasize their connection to the gravity side. 

^ Such a map also indicates that the black hole singularities are resolved in quantum gravity, 
which corresponds to the boundary gauge theory at finite A'' (the number color). See [11 1 for a 
discussion. 

^ The method was motivated from the approximation developed in [11| to compute boundary 
correlation functions. 
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Compared to many other approximations already existing in the hterature (see e.g. 
and [1T2| - ^ for an incomplete list), our method has the following features: 

• We introduce scaling limits which amount to a (generalized) geometric optics approx- 
imation. In the limits there is a natural connection between quasi-normal poles and 
complex geodesies of the black hole geometry as in 

• The quasinormal frequencies are expressed elegantly in a quantization formula of 
Bohr-Sommerfeld type. We believe the formula should also apply to perturbations of 
a vector field or metric perturbations in a generic black hole geometry (e.g. charged) 
in AdS, once one has written the equation of motion of a given field in terms of 
a Schrodinger equation. Similar quantization formulas have been derived before in 



the case of asymptotically fiat black holes starting with [|T2[ and for AdS black holes 
in ||2^ . Compared to these earlier work, our scaling limits introduce a small parameter 
and thus give the precise regimes to which the quantization formula can be reliably 
applied. This allows us to determine the location of quasinormal modes with very high 



frequency /overtone number which were not safely determined in previous work ||26 



Our method can also be applied to determine the spectrum of quasinormal modes for 
any complex frequency value, including quasinormal pole lines starting far from the 
real axis. 

• The method applies well to quasinormal perturbations with large angular momentum 
(or momentum) which are hard to address using other methods. 

Applying our Bohr-Sommerfeld type formula to an AdS black hole geometry dual to 
a CFTd on IR^''^"^ or IR x S'^~^, we find the the following results among others: (we will 
quote the results for d = 4, see main text for other dimensions) 

1. For a scalar operator O of dimensionB A ^ 1 in a CFT4 on IR^'"^, for momentum 
p/T <^ A, the quasinormal poles of the corresponding bulk scalar field are given by 

^ (±1 -z)7rT(A + 2n- 1) +0 ( ^ ) , n = 0, 1,2, ■■■ (1.1) 



^A, 

where T is the temperature. Equation ( |1 . 1| ) implies when the thermal system is 
perturbed by the operator (9, the relaxation back to equilibrium happens very fast, 
given by a time scale 



This corresponds to a bulk scalar field of mass m which is related to A by A = ^ + 



' ^ + m?R'^, where d is the spacetime dimension of the boundary theory and R is the curvature 
radius of AdS. 
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Also note that the spacing of the pole is independent of A. Similar results apply to 
a CFT on a sphere for angular momentum not too large. For k = ^ 1 we 

instead find that 

w^p±e^*(7rT)(4A;)"^ (^A -2 + 76(271 + 1)), n = 0, 1, ■ ■ ■ (1.3) 

Compared to ( p..2|) this implies that at large k, the relaxation time is enhanced by a 
factor k ^ . Here the order of limit is important in deriving (|1.3| ) we take A ^ oo limit 
first and then k large. 

2. For any A (not necessarily large), in the limit p ^ ttT, the quasinormal poles can be 
written in 1/p expansion as 

w = p±0.344e^'?(7rT)^p-^A-l + 2n)^ n = 0, 1, ■ • • . (1.4) 



Equations ( |1.3| ) and ( p..4|) have qualitatively similar behavior in their overlapping 
region of applicability, but not exactly the same, implying the two limits of taking A 
and p large do not commute. Note that for (|1.4|) the spacings of the poles are not 
homogeneous. For p ^ (A — l)^7rT, the poles in ( |1.4| ) have an imaginary part which 
is much smaller than the temperature scale, but the imaginary parts are always of 
the same order as the spacings between the poles. Thus in the spectral function we 
expect not to see individual peaks for each pole in (|1.4| ) , rather a broad peak close to 
the light cone co = p refiecting the effects of a large number of poles. In this sense 
( p..4|) are more like a branch cut than quasi-particle poles. Similar remarks also apply 

to(|iD. 

3. New interesting behavior arises for a theory on S'^~^ in the large angular momentum 
regime. Here one finds (as already noticed by [^) that for an operator of any di- 



mension, when the angular momentum / of a quasinormal mode is big enough, there 
exists a large number of quasinormal poles with a small imaginary part, proportional 
to e"'*- "-*. The spacing of the poles remains finite in the large / limit. Thus we find a 
large number of long-lived quasi-particles with a large angular momentum!!. The exis- 
tence of a large number of quasinormal modes with a very small imaginary part follows 
from the following simple fact of the bulk geometry. For a black hole of a spherical 
horizon, there exists a stable circular orbit for a particle moving in its geometry when 



Long-lived quasi-particles with a large momentum were also found in a CFT on ]R^'^ dual 



to a Gauss-Bonnet gravity in [28,29|. 
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the angular momentum of the particle is big enough, since the centrifugal force pre- 
vents the particle from falling into the black hole. While classically the orbit is stable, 
in the quantum theory the particle can nevertheless tunnel into the black hole, giving 
rise to a small imaginary part. From the boundary theory point of view, however, the 
existence of such long-lived quasi-particles for a theory on S'^~^ is curious. We discuss 
their possible physical origins at end of section 6. 

The organization of the paper is as follows. In section 2 we give a brief review of 
the geometry of a Schwarzschild black hole in AdS. In section 3 we rewrite the equation 
of motion in terms of a Schrodinger equation and discuss the condition for determining 
the quasinormal modes. In Section 4 we introduce a WKB approximation which can 
be considered as a generalized geometric optics approximation and present a quantization 
formula of Bohr-Sommerfeld type for finding the quasinormal modes. We apply the formula 
in section 5 to a theory on IR'^~^ and in section 6 to theories on S'^~^. The detailed 
derivation of the Bohr-Sommerfeld type formula is given in Appendix A. 



2. Black hole geometry 

In this section we briefly review the classical geometry of a Schwarzschild black hole 
in an AdS^+i {d > 3) spacetime. The metric can be written as: 

ds^ = -f{r)de + f{r)-^dr^ + r^dQl_^ (2.1) 

with 

/W = -' + l-^' (2-2) 
where ^ is proportional to the mass of the black hole, and d^'^_'^ denotes the metric on 
a unit {d — l)-sphere. We have set the curvature radius of AdS to be unity, as we will do 
for the rest of the paper. As r — oo, the metric goes over to that of global AdS. The fully 
extended black hole spacetime has two disconnected time-like boundaries each of topology 
S'^~^ X IR (see Fig. 1 for the Penrose diagram). The event horizon radius tq is given by 
the unique positive root of the equation 

+ 1 - = 

and the inverse Hawking temperature is given by 

When |U > 2, (|2.1|) dominates the thermal ensemble and describes a boundary CFT^ on 
S'^~^ X IR at a temperature given by 1/(3 [||. 
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Fig. 2: Penrose diagram for the AdS black hole. A null geodesic going from the 
boundary to the singularity is indicated in the figure. 



There also exists a black brane metric 

ds'^ = -fdt^ + jdr^ + r^dx^ (2.4) 

where dxf denotes the metric for a flat {d — l)-dimensional Euclidean space and 

f = r'-^- (2-5) 

The inverse Hawking temperature is 



The boundary manifold now consists of two copies of IR ' . The metric (|2.4| ) describes a 
boundary CFT^ on 1R^''^~^ at a flnite temperature given by (|2.6|). Note that for a CFT^ 
on IR^''^"^, since the temperature is the only scale, it can be normalized at any given value. 
We choose ( |2.6| ) for convenience. The boundary coordinates (t, Xi) are dimensionless with 



basic unit give by Note that (U) can be obtained from ( |2.1| ) by taking the following 
scaling limit 

ro ^ cx), r ^ ror, t t/ro, rldfl^-i = d^l ■ (2-7) 
This is the high temperature limit of the boundary CFT on S'^~^ . Equivalently we can 
flx the temperature scale as in ( p.6|) , then effectively the size of the sphere goes to infinity 
and we get a theory on IR'^"^. 

To describe the black hole geometry it is often convenient to use the tortoise coordi- 
nate: 

dr 1 

where are zeros of / with tq being the horizon. The region outside the horizon (region I) 
corresponds to 2 G (0, +oo). At the boundary r ^ oo we have z ~ ^ — > 0. At the horizon 
r ^ ro we have z ^ —j- log(r — tq) +oo. 
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3. Propagation of a scalar field in an ArfS'd+iSchwarzschild black hole 

Consider a scalar field! with quadratic action 

S = -\j dr(fx [id(p)^ + m^^] (3.1) 



in the background of (|2.1|) or (|2.4| ). According to the standard AdS/CFT dictionary, the 



dimension A of the boundary theory operator O dual to (p is given by 



A=- + z., ry=^ni^ + -. (3.2) 

For (U) let 

(P = e-'^^e'P-^r-'^'i(j{uJ,p;r), (3.3) 

the Laplace equation for (p can then be written in terms of the tortoise coordinate z ( |2.8| ) 

as 

{-d^^+Vpiz)-u^)i, = 0, (3.4) 
where Vp is an implicit function of z (below = p- p) 



Vp{z) = f{r) 



rf.2 ^ 4r^ 



(3.5) 



For ( p.l|) one replaces the plane wave in the x directions in ( |3.3| ) by spherical harmonics 
on S'^-^ and get Q with now the potential given by 

where / is the angular momentum on S'^~^ . 

As discussed below (|2.8| ), the region outside the horizon corresponds to 2 G (0, +00) 

with 2 = at the boundary and z +00 at the horizon. Both ( p.5| ) and (|3.6| ) behave 
near the boundary as 

Vp ^ ^, z^O, (3.7) 

z^ 

and near the horizon 

Vp oc e"^^ ^0, 2 ^ +00 . (3.8) 

The fact that for Re^; ^ 1, r is a one-to one periodic function of z with a period implies 
that Vp[z^ can be expanded for large real z as 

Vp{z) = Y,cir.e-'i^' . (3.9) 

n=l 

This property will be important in our discussion below. 



^ Since the background Ricci scalar is a constant, the term below should be considered as 
the sum of the standard mass term and the coupling to the background curvature. 



7 



V(z) 



v(z; 



Fig. 3: Plots of the potential (|^) for A; = 6, v = 2 (left) and with 
^ = 6, /i= ^, i/ = 2 (right). 



Note that the potential ( p.5| ) is positive definite and monotonic in z G (0, +oo) for 
z^^ > i and any > (see fig. 0). The potential ( p.6| ) is monotonic as in fig. ^ for / 
smaller than a critical value /c, while for / > the potential develops a well as in fig. ^ 
/c can be found by solving V'{r) = V"{r) = 0. Its explicit value is rather complicated and 
we do not give it here. An implicit expression in the large u limit will be given in sec 7.2. 
The potential well refiects the fact that when the angular momentum is sufficiently large 
there exist stable orbits for a particle moving outside the horizon. 



In our discussions below, we will use the notation appropriate for (|2.4| ). The discus- 
sions apply to both ( |2.4| ) and ( p.l| ) unless mentioned explicitly. 

For any given real uj, the Schrodinger equation ( ^.4| ) has a unique normalizable mode 
9w,p, which we will normalize near the boundary as 



0. 



(3.10) 



The asymptotic behavior of gujp{z) near the horizon {z — > oo) can be written as 



z ^ oo 



(3.11) 



which defines the function f{uj,p). Since the boundary condition (|3.10| ) is real and inde- 
pendent of i^iPi for general complex values of Wjp, we have 



r{u,p) = f{-u\p*) 



(3.12) 



When plugged back into ( |3.3| ) the first term in ( |3.11| ) describes a plane wave falling into 
the horizon at t ^ +oo, while the second term describes a wave coming out of the horizon 
at t ^ — oo. 
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At a given p, quasinormal frequencies for (p are those (complex) values of uj for which 
the normalizable solution is in-falling at the horizon, i.e. the solution behaves as e^^^ for 
z ^ oo. From (|3.11|) , the quasinormal frequencies iVq{p) thus satisfy 



/K,P) = 0. (3.13) 

One can show that if the potential ( |3.5D does not admit bound states (which is true for real 
momenta) the only zeroes of f{uj,p) are in the lower half w-plane. Also for real p it follows 
from ( p.l2| ) that if ujq is a quasi- normal frequency, so is — c<;*, i.e. they are symmetric with 
respect to the lower imaginary axis in the complex cu-plane. 

Using the prescription for computing the retarded Green function in AdS/CFT ^ the 
quasi-normal frequencies are precisely the poles of the momentum-space retarded Green 
function Gr{u!,p) for the operator dual to (j) in the boundary gauge theory. In particular, 
those whose imaginary part is much smaller then the real part correspond to quasi-particles 
in the boundary theory. 



4. Quasi-normal modes from WKB approximation 

We now develop an approximate method to determine the function gujp{z), from which 
we obtain the quasi-normal modes which are zeros of /(w, p) defined in (|3.11|) . This method 
is applicable to values of the momentum p of the same order as u) in contrast to other 
approximations considered in the literature [1T6| , [T7[] which require uj ^ p. We first introduce 



the approximation and then summarize our results for the quasinormal frequencies, leaving 
the details of their derivation to Appendix A. For simplicity of notations, below we will 
use the notation appropriate for the flat case ( p.5| ) . The discussion applies without change 



to the sphere case (|3.6| ). 



4.1. Large mass limit 

Consider the following large u limit with u and k fixed 

UJ = uu, p = fk, z/ ^ 1 , (4.1) 

i.e. we take the mass m of to be large and "measure" the frequency uj and momentum p 
in units of m. This is analogous to a geometric optics limit in which we take uj and k large 
but keep the velocities fixed. In this limit, one can solve ( p.4|) approximately using the 
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standard WKB method with l/u playing the role of U. More explicitly, writing ^ = e^^ 
equation ( 3.4|) becomes 



with 



and 



-{d,Sf - -dlS + V{z) + \q{z) = 



V{z) = / 1 + -2 



(4.2) 
(4.3) 
(4.4) 



V(z) 








Fig. 4: Schematic plot of the potential (| 
the horizon at z = cx). 



for k & R. The boundary is at 2; = 



We first restrict to real /c, for which case the potential ( [4.3| ) (see fig. ^ is a monoton- 
ically decreasing function for z G (0, +cxd), and consider positive energy scattering sates 
with u > 0. Solving ([4.2| ) order by order in l/u, we find that in the classically forbidden 
region (i.e. for z < Zciia. fig. ^ the exponentially decreasing solution to ( |3.4| ) can be written 
as 

, ... 4 „ , . 

(4.5) 

with0 

Z{z) = / dz' k{z'), 

J 2c 

(4.6) 





A 


e 

y/K{z) 


Z{z) = 


/ dz' 

J Zc 


k{z) = 




A = 


lim e"^^ 



^ The branch cuts for k{z) on the complex z-plane are chosen so that they do not intersect the 
integration contour in Z. 
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Zc = z{rc) in the lower integration limit of ( [4 .61) is the turning point with given by the 
real positive root of the equation 

V{r) = f{r) [l + ^^=u' . (4.7) 

For /c^, v? > 0, equation { \i.7\ ) has a unique positive root Tc > 1. 2{z) satisfies the equation 

1 k'^ 1 

7^'' - ^ + 7^' = 1 (4-8) 

with Z'{zc) = 0. Note that to the order indicated in ( [4.5| ), the term in (|4.2|) proportional to 
Q(z) does not contribute. As discussed in detail in [jll| in a classically forbidden (allowed) 



region equation (|4.8| ) is simply the geodesic equation for a space-like (time-like) geodesic. 

The expression for gl^^^ in the classically allowed region of the potential ( ^TS] ) (i.e. 
for Zc < z < oo) follows from the standard connection formula 

= cos (z.M^ - 4 ) (1 + O {y-')) (4.9) 

with 



l^z'piz'), p(z) = Ju^-f(^l + ^'^ (4.10) 



Near the event horizon (|4.9|) has the form 

A 



glf^ = e'^'+'^- (1 + O (z/-i)) + c.c. (4.11) 



Higher order 1/z/ corrections in (|4.5[ ) and (|4.9| ) can also be obtained from ( |4.2p us- 
ing the standard WKB procedure. In particular, the term proportional to Q{z) will be 
important at order . 

4-2. Large momentum limit 

For a scalar field with a mass v ~ 0(1), it is possible to use WKB methods to study 
( p.4| ) (or ( |3.6| )) in the large momentum (or angular momentum) limit p ^ oo with u fixed. 
Define 

uj = pu (4.12) 
then for p ^ oo equation ( 3 -41) becomes for tp = e^^: 

-af ^ + Viz) + i 
p p^ 



-{d,Sf - -dlS + V{z) + -Q{z) = (4.13) 
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with 



V{z) = ^ (4.14) 



and 

ew = /(^ + ''^-i) ■ 

When z (i.e. r oo) the potential ( [4.14|) is constant and equal to 1, which is different 
from (|3.7|). As a resiilt the WKB solution near the boundary behaves like: 



which is not the correct behavior (|3.10|) . This can be fixed as follows. One can solve the 



full equation (|3.4|) for g{z) near z —>■ and then match it with the WKB solution ( [4.16|) . 
For u > 1, solving (|3.4D near the boundary z ^ where the potential V{z) ~ 
+ z~^(z/^ — i), we find that 

g^p{z) oc z^J,{p{u^ - . (4.17) 

In the large p limit, for any z ^ 0, we can expand ( |4.17|) for large pz, giving 

g^p{z) ~ cos (^piu^ - l)h - ^ (J + 0) (4-18) 
which can be matched with the WKB solution and fixes g^^^^^ to be 

g^jJ-'Hz) ~ cos (pW - ^(1 + 2Z.)) , (4.19) 

with 

W = f dz'p{z'), p{z) = ^v? - f{r)r-^ . (4.20) 
Jo 

Thus for u > 1, we get parallel result with ( [4.9| )-( |17TO[ ) with z = as the turning point in 



place of Zc- Also note there is an additional ^ term inside the argument for cos in ( ^.19| ) 
which differs from ( [4.9|) . 

When u < 1 there is now a turning point at z > 0; apart from possible exponentially 
subdominant termsl the result is exactly the same as ( |4.9| ) with p{z) = ^Jv? — f{r)r~'^ 
and p in place of u as the large parameter. 



^ For u with a small imaginary part these can be obtained by a more careful WKB analysis as 
outlined in Appendix A. 
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4-3. A Bohr-Sommerfeld quantization formula for quasinormal frequencies 

The positions of the quasinormal modes in the large u or p limit can be obtained by 
generalizing the WKB analysis of the previous subsections for complex values of u. A very 
general but somewhat involved procedure is described in appendix A. Here we outline the 
main ideas and give the final results. 

The quasinormal modes are zeroes of the coefficient f{uj,p) of the second term in 
In the WKB expansion, as is clear from equation f{uj,p) can have zeros 

only if certain exponentially subdominant term becomes of the same order as the dominant 
term ( |4.11| ) , and the two cancel each other. The values of u at which certain subdominant 
terms become of the same order as the leading term form dimension one lines in the 
complex-w plane, the so-called anti-Stokes lines of the expansion. The quasinormal modes 
are then special points on the anti-Stokes lines at which the two contributions exactly 
cancel each other. 

For real u the WKB approximation of the function Qujpiz) for large z is given by ( [4.9| ). 
To find quasinormal modes we first analytically continue ( [4.9| ) to lower complex plane. We 
start with large positive values of u and after giving u a small negative imaginary part, we 
can rewrite ([4.9|) asi 

where e is an IR cutoff near the boundary and we have expressed the normalization factor 
A in ( |4.9D explicitly using ( [4.61 ). In the first term the integral in the exponential of A in 
( [4.6|) and that for W in ( ^4.9| ) combine into a single contour integral from e to z, which 
does not depend on the turning point. In contrast the second term depends explicitly on 
the turning point Zc{u), since we have to evaluate two integrals one from e to the turning 
point Zc{u) and a second one from Zc{u) to z. Near the horizon z oo, the first term 
in ( [1.21| ) behaves as e*'^^^ while the second as e"**^"^, which can be identified respectively 
with the two terms in ( |3.11| ). Note that had we given a positive imaginary part to u the 
role of the two terms would be interchanged, i.e. it would be the term behaving as e^'^^^ 
that depends on the turning point Zc{u)- 

As one varies u in the lower half plane, there are special values Ub, at which Zc{u) 
merges with some other turning point of equation ([4.7| ), which we denote as zt{u). More 

^ For a square root function like the one in k{z), we take the branch cut to be along the 
negative real axis. 
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explicitly, at Ub, we have zdub) = zt{ui,) = Zb and V'{zb) = 0. In the case that zt{u) is an 
active turning point, it gives a subdominant contribution (except on the anti-Stokes line 



specified below) to gl^p of exactly the same form as the second term ( |4.21| ) with Zc replaced 



by zt- There is no correction to the first term in ( [4.21|) , since it does not dependent on 



the turning point. Thus in the large z region, the ratio of the contributions from Zc and 
Zt to f{u!,p) in ( |3.11|) is given by 



The two contributions are of the same order when Z{zc, zt) k{z') is pure imagi- 

nary. This condition determines an anti-Stokes line in the complex u plane passing through 
Ub, when the ratio becomes unity. More careful WKB analysis shows that the anti-Stokes 
line is in fact a half line since along the other half line zt ceases to be an active turning 
point and needs not to be considered. Which half line to be chosen can often be easily 
determined on physical ground without detailed analysis. For real spatial momentum, the 
anti-Stokes line should always emanate from Ub in the direction away from the real w-axis, 
since for real u there is no exchange of dominance and thus the anti-Stokes line should not 
intersect the real w-axis. Note that Ub is also a branch point for Zc{u) and in our discus- 
sion above we have assumed that Zc{u) is always the dominant turning point outside the 
anti-Stokes line, which is equivalent to choosing its branch cut to be along the anti-Stokes 
line. 

f{uj,p) in (|3.11|) has zeroes when q'^^z{zc,zt) _ _]^^ which leads to the condition 

2uZ{zc,zt) =2u dz'K{z') = ni{l + 2n), n = 0, 1, ■ ■ • . (4.23) 



Note that only half of the integers are chosen in ( |4.23| ), refiecting the half line nature of 



the anti-Stokes line. Equation ( [4.23| ) determines a line of quasinormal poles which starts 



at Ub- To find all quasinormal poles we can repeat the procedure for all points where Zciu) 
merges with some other turning point. 

To summarize, we now give the full procedure for finding the quasinormal modes 
using (|;2|): 

1. Find the turning point Zc{u) from ( ^.7| ) for large positive u and analytically continue 



it to the lower half w-plane. Note that while ( [4.7| ) have multiple solutions Zciu) can 
be uniquely determined for a given large positive u. We will call Zciu) the physical 
turning point. 
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2. Find the roots Zf, to equation dzV{z) = and from the turning point equation 
find the corresponding Ub = V{zh). At such a Ub two turning points of {\i.7\ ) merge 
together. However not all of them are relevant for finding quasinormal modes. One 
should only consider those Ub while lie in the lower half plane and correspond to the 
merging of the physical turning Zc{u) with some other turning point. We will call such 
a Ub, physical Ub- 

3. For each physical Ub, there is a line of quasinormal poles determined by ( |4.23|) , where 
zt{u) is the other turning point which merges with Zc at Ub- The direction of the line 
should point away from real w-axis. 

Equation ( [4.23| ) can be greatly simplified for the lower-lying quasinormal modes close 
to Ub- Let u = Ub + X with x ~ 0{u~^). Near Zb = zriub) = Zc{ub) we can approximate 
the potential V{z) as 



V{z)^V{zb) + ldyizb){z-Zbf. 



We then find that 



..,.«..±a, «-2./5p^ (4.24) 



and after integrating (f4.23| ) 



Ziz^,zr) - ±'l^Jm^ = (4.25) 



with 



2V 



(4.26) 



Equation (|4.23|) then leads to the position Un of the quasinormal modes 



1 S 

Un=Ub + in+ -)-, n = 0,l,--- (4.27) 

In terms of w = uu, the modes are uniformly spaced near Ub — i^Ub 

1 

Un = ujb + {n + -)d, n = 0,l,--- (4.28) 

with a spacing given by0 ( f4.26| ) which is independent of z^. The branch of the square root 
in (|4.26| ) should be chosen so that 6 points away from the real axis at Ub- Note that for 



The formula below does not apply to the non generic case d1V{zh) = which is reahzed if 
more than two turning points merge together or if V{zb) = for finite Zb- 
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modes which are of order 0{u) from ui, in the w-plane, ( ^4.28| ) is no longer vahd and one 
needs to use ( [4.23|) . 



For large p and uj but finite u we can use the limit ( [4.12| ) to obtain expressions for 
the quasinormal modes. We will start from large real u and analytically continue in the 
lower half plane. As described earlier the boundary z = should be considered as the 
turning point for u > 1. At zi = 1 this exchanges dominance with (zt) present for u < 1. 
An anti-Stokes line starts from u = 1. On this line we find poles for: 

2pZ{0,zt) =in{l + u + 2n), n = 0, 1, ■ ■ • (4.29) 



where the difference with respect (|4.23|) comes from the extra ^ in (|4.19| ) with respect to 



5. Applications to a CFT on R^'''"^ 

We now apply the results of previous section to find the quasinormal frequencies of 
scalar perturbations of an AdS black brane with metric ( |2.4| ), which correspond to poles 
of retarded Green functions for a CFT on IR^''^~^ at finite temperature. 

Let us first consider k — 0, in which case 

V{z)=firiz))=r'-^. (5.1) 



The turning point equation ( [4.7| ) has coincident solutions for those values of r such that 
V'{r) = that is 

2.+ ^ = 0. (5.2) 

Among d solutions of ( |5.2|) only the following correspond to physical ones, i.e. where the 
physical turning point Zc{u) merges with some other turning point, 



r = — - — e 



The corresponding line of quasinormal modes starts at 



ui, = ±Cde^'^, cd=(^—) \ -^. (5.3) 



d-2 
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The spacing of modes near uj, can also be easily computed from (^4.26| ) 

6 = ±Vdcde^'^, (5.4) 

where each sign of (5^) should be paired with that of (5^). Expanding ( 4.23|) for large u 
we can find the mode spacing as |c<;| — oo: 

5 = ±dsin^e"^'^ (5.5) 

Except for (i = 4 (i.e. AdSs), ( ^.4| ) is different from the spacings (|5.5|) at large u 
although they do point to the same directions. This is expected as for k = since from 
( ^.1|) that the phase of V{z) is constant for arg{z) = and therefore the mode line is 
straight. As we will see this changes for /c 7^ 0. 

The explicit analytic expressions for the branch points and spacings of the quasinormal 
modes for k are rather complicated since they involve roots of high order polynomial 
equations. We will only point out some important features for k large. 

For a real k, the location of the quasinormal modes is similar to that for k = 0. There 
are four lines of modes. In particular, as /c^ +00, the branch points approach the real 
axis 

/ 2 

' " zn 2-n-i _ d-2 . 

e+d+2 k d+2 I (5.6) 




and the mode spacings near the branch points are 

'd 



S ^ ±Vd + 2 



^2jri_ d-2 
e+d+2 k d+2 -f 



(5.7) 



One can check that the angle between the lines of quasinormal modes and the real u axis 
increases monotonically as /c^ increases. At large uj the effect of a finite k can be neglected 
and the mode spacing is again given by (|5.5| ) therefore the mode lines are no longer straight. 




Fig. 5: Numerical determination of one of the quasinormal mode lines in the 
complex u plane for k = 20, the straight lines tangent at 00 and at the starting 



point are obtained from (|5.5|) and (5.7) respectively. 
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In fig. ^ we plot a numerical determination of a mode line for k = 20 in d = 4 and 
the value of ( [4.23|) along this line as a function of the real part of u. 

At finite u, we can use the limit ( [4.12| ) to obtain expressions for the quasinormal 
modes at large p and w. As described there in this case z = (the boundary) should be 
considered as a turning point. The second turning point is 

rj^ = (1 - -u^) 2^ 

and merges with z = for u = 1. Equation ( |4.29| ) then gives the lower- lying quasinormal 
modes 

UJ = ±p±Ke'^^^ + + 2n)^, n = 0, 1, ■ ■ ■ (5.8) 

with K a numerical factor given by 

^ _ V^T[3/2 + l/d] 

r[i/d] 

Note that in ( p.8| ) the spacing between various modes increases with n for d > 2. In the 
limit w — > oo, one finds from ( [4.29| ) that the spacing between quasinormal modes tends to 
the limit 

5 = ±(isinQ) e^'^ 

which agrees with formula ( |5.5| ). 

Note that in deriving (|5.6| ) and (|5.7D , we first take the large u limit and then k large, 
while (|5.8| ) was obtained by taking p ^ oo with u fixed. Taking n = and extrapolating 
( ^.8|) to the regime p = ku with u large, we find it gives an expression which agrees with 

(after multiplying ( |5^ ) by u) up to the factor of K. Thus we see the two limits do 
not quite commute. 

We now summarize our results by writing down various expression for d = 4 explicitly. 
We will only write down the sequence on the lower right quadrant. The sequence on lower 
left quadrant can be obtained by refiection with respect to the imaginary w-axis. From 
)- (|5.5|) , in the large u limit with k — {or small /c = ^ -C 1), 

{l-i)nT{u + 2n + l), n = 0, 1, 2, ■ ■ ■ (5.9) 



where we have restored the temperature dependence explicitly by using T = ^. 
differs from expressions in the literature [jl^ valid at large u which can be obtained by 



replacing n ^ n + ^ log(2). This is no surprise as the approximation used in |jT^ is not 
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valid for ~ a; ^ oo. The approximation of which uses a combination of large v and 



/c = also results in (|5.9|) . We expect (|5.9|) to be reliable for large values of uj as the regime 
of validity of [1^] is c<j ^ v^~^ . 

In the limit of large v and /c = ^ ^ 1 , (|5.6| ) and (|5.7| ) give the lower-lying quasinormal 
frequencies 

u^p + 2-h-^{7tT)(^-£^^ ' (^z/ + V6(2n+1)) , n = 0, 1, ■ ■ ■ (5.10) 

In the limit p ^ u with u fixed, we find from 

uj = p + Ke~'S- ('KT)^p-^ {l + iy + 2n)^, n = 0, 1, ■ ■ ■ (5.11) 

with K = 0.344. Note that ( ^.11|) also applies to massless fields in the bulk which corre- 
spond to = 2. 

For p ^ (A — l)^7rT, the poles in (|5.11|) have an imaginary part which is much smaller 
than the temperature scale. However, in ( ^.11|) the imaginary parts are of the same order 
as the spacings between the poles. Thus in the spectral function we expect not to see 
individual peaks for each pole, rather a broad peak close to the light cone uj = p reflecting 
the effects of a large number of poles. Note that the limit in going to zero temperature 
is not continuous, since for a CFT on IR"^ at flnite temperature, the temperature sets the 
scale for the system. At zero temperature the pole line (|5.11|) is replaced by a branch cut 
starting aX ui = p. 

The emergence of the time scale p^/^ as implied by the imaginary part in (|5.11|) can 
also be seen by looking at time-like geodesies propagating in the bulk. For a time like 
geodesic we have 

— = ^ fir) = - — (5 12) 

dr f(r)y/E^ - f{r){l+pyr^)' 

where E and p are the conserved quantities along the geodesies corresponding to momen- 
tum in t and in the spatial directions. For large p, equation ( |5.12| ) has the following scaling 
limit 

E'^-p^r^P^, r P^ , t P^ . (5.13) 

Thus such a geodesic approach the boundary more and more as p is increased and the time 
the geodesic spends in the region where r is 0{p^) is also of order p^ . During this time 
the geodesic propagates close to the boundary light cone. 
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To conclude this section, let us briefly comment on the quasinormal modes for pure 
imaginary /c's. Imaginary /c's are of interest as they play a role for probing the geometry 
inside the horizon and near the singularity using Yang-Mills theory |TI|]. For /c^ = —q^ < 0, 
there is an additional line of quasinormal modes along the imaginary u axis with branch 
points located at uq = —iEc- For <^1, one flnds that 

, _ d-2 

Ec^\-A^r^] (5.14) 



and the spacing between modes near the branch point is given by: 



d 



Ec decrease monotonically to zero as increases to 1. For ~ 1 the branch point reaches 
M = and the two merging turning points tend toward 2; = +00. The mode spacing can 

d.V{z) I 



still be found using (|4.26| ) because lim^^^o y^^) Uo is well defined. This results in a 



spacing 5 between successive modes 5 = id + 0{q — 1] 




Fig. 6: The potential ( |4.3| ) for —k'^ = > 1 admits bound states. In the z 
coordinate the boundary is at 2; = and the horizon at z = 00. 



For q^ > 1, this line of modes along the imaginary w-axis crosses the real axis. This 
is due to the fact that for q^ > 1, the Schrodinger problem (|3.4| ) admits bound states, as is 
clear from the shape of the potential plotted in fig. |^. The modes lying in the upper half 
plane are given by the bound states. In fig. 0, we plot schematically the configurations of 
the quasinormal mode lines for various values of /c^. 



6. Long-lived quasiparticles for a CFT^ on S'^ ^ 
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Fig. 7: Schematic structure of the mode hues for (a): = 0, (b): > 0, (c): 
— l<k^< 0, (d): fc^ < —1. At finite u, the modes are at discrete points along the 
fines. 



V(z) 





■^min ^max Z " 

Fig. 8: The left figure is a schematic plot for the potential V{z) for k > kc- The 
resulting quasinormal mode lines are shown in the diagram at the right. In this 
plot we only show the right half of the complex u-plane. The modes in the left half 
are obtained by reflection with respect to the imaginary axis. 



In this section we consider quasinormal modes for black hole metric and ( |2.2| ), 
which correspond to poles of retarded functions for a CFT on S'^~^ at finite temperature. 
In this case the quasinormal modes are labeled by angular momentum / on S'^~^ and the 
corresponding large u limit is now 

uj = i^u, {2l + d-2) =2iyk, > 1 . (6.1) 
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For small values of A; G IR, the potential V{z) ( [4 .31 ) is a monotonically decreasing function 
along the real z axis. The pattern of quasinormal modes is very similar to that of the 
black brane discussed in the last section and is qualitatively captured by fig. 0. Here we 
simply note that for d = 4 at k = (or /c ^ 1) one gets from ( f4.23| ) the following simple 
formulaO (for the sequence in the lower right plane) 



with B given by 



2tt 

= —{iy + 2n + l) n = 0,l,2,... (6.2) 



2tt 

B= , with lmB = (3 (6.3) 

ri — zro 



where ro is radius of the horizon and ri is defined b 



The other sequence of modes are obtained from the refiection of ( |6.2| ) across the imaginary 
axis. ( |6.2| ) is consistent with previous results in the overlapping region of validity. 



For larger k however there exists a critical value kc such that for k > kc the potential 
is no longer monotonic and looks like the left plot of fig. ^. At k = kc we have Zmin = Zmax 
and aX z = Zmin = Zmax both the first and second derivative of V{z) are zero. By using the 
relation ( |2.8| ) between z and r and the explicit expression of V{z) ( [4.3| ), kc can be found 
from the largest positive root of the coupled equations V'{z) = and V"{z) = 0. 

For example, for d = 4 one finds 

{kl - ^if - 21 = . (6.4) 

For large ( |6.4| ) give^l kc ~ 3-\/3a*- For generic dimensions d > 2, one finds that for 
large 

2 

kc ~ . (6.5) 



The form of the potential fig. ^ for k > kc implies that it is classically possible for a 
particle with sufficient angular 
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For other dimensions, the expression is more complicated. For d = 6 and d = 3 it is possible 



to evaluate ( 4.23 ) for = in terms of elliptical integrals. 

Note the mass of the black ^ can also be written as // = VQrf. 

which is pretty good already for /i > 2, required for the black hole solution to dominate the 
thermal ensemble. 
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momentum to be in a bounded orbit outside the horizon of the black hole as the 
centrifugal potential provides a barrier to its falling in the horizon. Quantum mechanically 
however there will be a nonzero probability for the particle to tunnel through the centrifugal 
potential barrier and be absorbed by the black hole. 

Applying ( [4.23| ) to the situation for k > kc (fig. ^) simply reduces to the standard 



Bohr-Sommerfeld quantization in the potential well (see also for the case of a massless 
perturbation in d = 3). We thus find resonances whose energies are given by uin — i^Un 
with Un determined by 

2u ' dz^/ul - V{z) = 27r (^n + , n = 0, 1, • • • (6.6) 

where w^^^ = V{zmin) < < V{zmax) = u^max ^nd Za, Zb are real solutions to V{z) = 
(see fig. P). The maximum energy for these quasi-stable states is Umax- For small n the 
first few energy levels will be close to V{zmin) and using ( [4.26|) we find that 

Un = Umin + ];\l ^2V{z' ^ ^ ^) ' ^ = 0, I,---, ul,^^ = V{Zrmn) ■ (6.7) 

These quasi-stable states can tunnel through the potential barrier between Zb and Zc (see 
fig. ^) to fall into the horizon, with a decay rate given by 

= exp (^-2u ^ ' dz^/V{z) - ul^ . (6.8) 



Note that ( |6.8|) is not captured by ( 4.23 ) since it is exponentially small in u. 



Equations (|6.8| ) and ( |6.7| ) thus imply the existence of poles with a very small imaginary 
part in boundary retarded Green functions. The quasinormal poles form lines in the 
complex w-plane extending on the real axis from Umin to Umax- At Umax the lines start 
deviating from the real axis as depicted in the right plot of fig. |[ For large \u\ ^ Umax 
the position pole line is qualitatively the same as the one for a monotonic potential. 

The large u limit is not essential to the existence of these resonances. For large 
angular momentum and finite u we can work in the limit 

uj = pu, {2l + d-2) = 2p, p>l (6.9) 

for which the WKB potential (|4.14|) is 



y(^) = 1 + 1-4 (6.10) 
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V{z) has a maximum at (see fig. 



1 



/ dfj, \ ""-^ ^ , , f d-2 \ , , 

rmax = \ ~^ j ' ^^^^ V{rmax) = ( ^ + ) (S-H) 



max 



For a given 1 < u < V{rmax), there are long hved quasiparticles for 

/■Zb 

2p dz^ui - V{z) = n{2n + 1 + z^), n G 0, 1, 

^0 



(6.12) 



with a decay rate 



= exp (-2P ^ ' dz ^/V{z) - ul^ (6.13) 



where Zb and Zc are the turning points before and after the maximum as labeled in fig. 
These long-lived quasiparticles extend from w = 1 to tt^ = V{zmax)- For w, ~ 1 



Ur,=p+{2n + l + u)\l^^=l + 2n + A, n = 0, 1, ■ • • . (6.14) 



Note that for equation ( |6.14D to apply, we need p{\/V{rmax) ~ 1) ^ For rmax large (i.e. 



jj, large), this condition can be written explicitly as 

/>A;U^~AT^ (6.15) 



where we have used that fx ^ for large fi. Note that ( |6.15| ) has the same scaling with 
fi as ( |6.5| ). 

Notice that ( |6.14|) is exactly the spectrum of a scalar conformal operator of dimension 
A with angular momentum / at zero temperature. That in the limit / — oo we recover the 
zero temperature spectrum is not surprising. We expect states of energies in ( |6.14|) not to 
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be excited in the thermal ensemble due to the exponential suppression in e ^' and that the 
degeneracy of these states only increases with I as a power. However, it is rather curious 



the power in T in the onset condition (|6.15|) and (|6.5| ) for the appearance of long-lived 



quasi-particles is much larger than 1, as one would naively expect. For example for d = 4, 
( |6.15|) gives / ^ AT^. Another interesting feature is that for a given I satisfying ( |6.15|) , 
the zero temperature spectrum ( |6.14[ ) only persists for a finite range of frequencies from 

2 

ujq = I + A to Umax ~ ^(1 + When uj > uJmax, finite temperature effects again set 

in. It would be interesting to understand these issues better and to investigate whether 
there are some physical processes in which such long-lived quasi-particles play a dominant 
role. 

To conclude this section, let us consider what happens to these resonances in the 
limit (|2.7| ), which describes the boundary theory on IR*^"^. In momentum space, the limit 
can be described as 

— = finite, — = finite, T ^ oo . (6.16) 



In the limit T — ^ oo, from equations (|6.15D and (|6.5|), the onset value of I for appearance 



2d 

of long-lived quasi-particles scales with T as T^-a. It then follows that the frequencies 

2d 

and angular momenta of the resonances scale with T at least as fast as T^-a > jp^ which 
is much faster than ( |6.16|) . Thus we conclude that these resonances disappear in the limit 
(p77|). Indeed, as we discussed earlier, the potential V{z) for (|2.5| ) is always monotonic and 
there is no kc- 
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Appendix A. WKB analysis of the quasi-normal modes 

In this section we describe how to extend the WKB approximation to account for 
subdominant contributions to g^jpiz) which are responsible for the appearance of the quasi- 
normal frequencies. We start by reviewing general techniques which are applicable to any 
WKB computation [PD|,PT[| . We then apply them to our specific case pinpointing some gen- 



eral features of the WKB result stemming from the structure of the potential ( |3.6| ) or ( |3.5| ) . 
Finally we work out a couple of specific examples. 
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A.l. Review of complex WKB 

By performing the rescaling ( |4.1D equation ( |3.4| ) in the large z/ limit reduces to: 

1 (f 



- V{z)+u^ ip{z) = 



(A.l) 



we want to find the form of g{z) defined in ( |3.1C1| ) for z ^ oo. From there comparing 
with ( |3.11| ) we will obtain approximate expressions for the quasinormal frequencies. 
Define Z{zo, z) = J^^ dzK{z, u) where k{z, u) = ^JV{z) — 




Fig. 10: Pattern of Stokes lines near a turning point 



Close to a generic point zq the lines of constant imaginary part for Z{zq,z) do not 
intersect. However whenever zq = zt with V{zt) = there are three of these lines0 
converging towards ^j- at a ^ angle; these are called Stokes lines. On each of these 
Z{zt, z) & R and so the exponential factor k{z, the WKB approximation 

to the solution is real and either decreasing or increasing along the line. Both K{z^u)~i 
and Z{zt, z) are multi- valued functions around zt and we will define them by introducing 
a branch cut extending from zt in region // in fig. [l^ . Suppose now Z{zt, z) < along 
line 1; then in the WKB approximation of the solution decreasing along 1 only the term 
w)~^e'^^*^^^'^) is present in regions / and //. We will find the WKB expansion of 
this solution in region /// by applying the principle of exponential dominance stating 

More than three lines is not a generic situation and is not considered 
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that crossing a Stokes line the coefficient of the dominant term in the expansion does not 
change. The term k(z, the dominant one on line 2 and so its coefficient 

is continuous and doesn't change while crossing it. However in region /// we can also 
have a contribution proportional to K{z,u)~'^e~^^''^'^'^^ which will be the dominant one 
along line 3. In order to find its precise coefficient we analytically continue the expansion 
we have in region // up to line 3 but in doing so we cross the branch cut and so we get 
-m(2,M)-^e-''^(^T'^). Then in region /// the asymptotic expansion is: 

K instead the branch cut was in region / we would have obtained in analogous way: 

Given a turning point zq its Stokes lines cannot intersect unless in another turning point 
which is not a generic situation and is excluded in what follows 111, and therefore divide 
the complex z plane into three regions. Two cases are possible: 

a) The origin and +oo are in the same region then zt is called inactive. 

b) They are in different regions and zt is an active turning point. 

For each active turning point imagine shading the regions in which the origin and +oo 
are; the intersection of all the shaded regions will be called active region in the following. 
Starting from the region containing the origin we can order the regions which are active 
by adjacency up to the one containing +oo. 

For every turning point there is a region which doesn't contain neither or +oo. It 
is therefore possible to choose the branch cuts defining \JV{z) — and Z in such a way 
that they do not cross into the active region. 

Then in the active region we can globally define two wave-forms exp(±z^2^(2;o, 2;)). 
For each wave-form the Stokes lines where it is dominant or subdominant are fixed by the 
sequence of active turning points as shown in figure. 

these non generic cases are the ones that divide the parameter space in regions with topo- 
logically distinct patterns of Stokes lines 
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Fig. 11: schematic representation of the active region. The red dashed Unes are 
branch cuts while the arrows are in the direction of decreasing real part for Z 

The WKB approximation to the desired solution is then given by choosing the appro- 
priate coefficients of the two wave-forms in each region. When we cross one of the Stokes 
lines the coefficient in front of the wave-form which is dominant along that line does not 
change while the coefficient of the other wave changes according to the appropriate con- 
nection formula. Two cases are possible and are shown in figure (the arrows represent 
decreasing Zi^zt.zy): 



Suppose we start with the wave Aexp(zyZ(zo, -z)) + -Bexp(— z^Z(zo, 2)) in region / 
then in the first case the first term is dominant while crossing the Stokes line and so its 
coefficient doesn't change. By applying the cormection formula we get in region //: 

Aexp(zyZ(zo, 2)) + -Bexp(— z^Z(zo, ^)) + ^Aexp(^^Z(^o, zt) — vZ^zt, z)) 
while in the second case the dominant term is the second one and we obtain: 

Aexp{uZ{zo, z)) — iBex.p{—uZ{zo, zt) + i'Z^zti z)) + -Bexp(— //^(^Oi z)) 

A. 2. General remarks to the Black hole case 

The application of the method just described to the specific case of the black hole 
background is in principle straightforward once the pattern of stokes lines corresponding to 
the values of k and u of interest is determined. Unfortunately the range of possibilities is 
quite extended and we will content ourselves to determine what are the turning points that 
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give the dominant contribution to the result in various regimes, and determine the position 
of quasi-normal modes. In this subsection we will comment on some general features of 
the WKB computation which arise from the exponential behavior of V{z) for z —>■ oo. In 
particular this implies that the boundary z — and the horizon z — oo are separated by 
an infinite number of Stokes lines. 

Being V{z{r)) a single valued function of r it follows from the discussion of the 
function r{z) in Appendix B that V{z) will have branch points at those points z for which 
r{z) = 0. In the following we will use the determination for r{z) described in Appendix B 
which has the following properties: 

1 Re{z) +00 corresponds to r{z) approaching the horizon, z = corresponds to the 
boundary. 

2 The only branch points present for Re{z) > are located at z = j + i {n + 
where tt, G Z and + i(3) is the location of the singularity in the complex z plane. 
For any branch point the branch cut extends on a line of constant imaginary part for 
z towards Re{z) = — oo. 

3 This determination is such that V{z + i^) = V{z) 




We will describe the consequences of the periodic structure of V{z) with the simple 
case arising for u close to the origin so that there is only one active turning point T (and 
its periodic images), which is the one considered in (|4.7| ) for G R'^. The pattern of 



Stokes lines is represented in fig. |T2| for — f < Arg{u) < 0. 
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Requiring Z to be continuous in the active region fixes it up to a sign. In the foUowing 
we wiU conventionaUy choose the determination of Z in the active region in the following 
way: For any two points Z\, Z2 = Zi + the quantity Z{zi,Z2) must be the same due 
to the periodicity of the function V{z). By computing it for Re{zi) +00 we obtain: 
Z{zi, Z2) = ±iu JJ^ dz. We will choose the sign of Z such that Z{zi, Z2) = — f 

For z in zone / in figure, the WKB wave-function corresponding to gujp{z) defined in 
is: 

9lo,p{^) ~ AK(z,'u)"2exp(z/Z(2T,2)) 

(A.2) 

A = lim exp(z^Z(e, 2;T))e'' = exp(z^Z*(0, 2;t)) 

with Z*{0,z) = lim,^o(2(e,^) +z^log(e)) 

Applying the rules given above we get for z G //, ///, IV. 



gcop{z) r^AK{z, u) 2 [exp{uZ{zT, z)) + iexp{-iyZ{zT, z))) 
Qujpiz) ^Ak{z^ u) 2 (exp(z/'2(2r, z)) + zexp 

9lop{z) r^AK,{z, u)~^ (exp(z^Z(2T, z)) + zexp(-z^Z(2T, z))[l + exp(2z^Z(2T, -ZtJ) + exp(2z/Z(2T, ^Ta))]) 
For z —>■ +00 we obtain: 

1 

1 - eyip{2uZ{zT, zti)) _ 



9ujpiz) ~ Ak{z, u) 2 I exp{h'Z{zT: z)) + iexp{-h'Z{zT: z)) 



we can use the fact that as zt-^ = zt + the quantity exp(2z^Z(zT, -^Ti)) = e to get 



' 2 



_i / \ 
g^p{z) Ak{z,u) ^ \^exp{uZ{zT,z)) + i^^-^^^jj^^-^j^exp{-uZ{zT,z))j (A.3) 

uuf3/2 

The factor 2sinh(i/u/j/2) '^hich is required by the analysis in Appendix B arises naturally 
due to the periodic nature of the potential. 

For z —>■ +00 we have Z{zt, z) ~ iuz and the following holdS: 



It would seem that this expression is inconsistent with the general form of gujp{z) described 
in Appendix B as the coefficients in front of e*"^"^ and g^*"^"^ are not related by n ^ —u. However 
the form of the WKB approximation to gu;p{z) does not have to be continuous in u. In fact ( lAlD 
was found for Im{u) < while for Im{u) > the following expansion is valid 

g^p{z) Ak{z,u) 2 \ eyi.p{-vZ{zT,z)) - I ^^^^^^^^^^^^ eyiY>{uZ{zT,z)) \ . (A.4) 
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^ — e--^+-^ + z-^-^^^e— — ^ . (A.5) 
ipu \ 2sinn(z^u/5/2) / 

where 

7= lim {Z{zt, z) — ifuz) (A. 6) 

z — >+oo 

Following our general discussion the quasinormal modes for Im{u) < show as zeroes 
in the coefficient multiplying e~^^'^^ in ( |A.5|) . We conclude that there are no quasinormal 
modes for small complex uj = uu. 

Due to the periodic nature of the potential however even in this simple case an infinite 
number of subdominant contributions to the WKB expression have to be considered for 
z ^ oo. This same complication will arise when other turning points enter the active 
region. 



A. 3. The k = case for large u 

Increasing u the pattern of Stokes lines changes as a second turning point becomes 
active. The situation is exemplified by the figures below representing the Polya vector 
field for the functionS K{r,u)f{r)~^ for = 4, d = 8 and various values of u. The 
singularity is represented by a purple dot at r = while the horizon is a green dot. In 
the first figure which is representative of the situation for small u there is only one active 
turning point represented in red. As the norm of u is increased a second turning point (in 
blue) enters in the active region (figure 2). For u close to the real axis the first turning 
point is the dominant one (figure 3), however as the phase of u becomes more negative 
the two turning point exchange dominance (figure 4) and finally it is the second turning 
point which is dominant for even more negative phase of u (figure 5). For |w| ^ oo the red 
turning point which we will denote as T approaches the boundary — > while the blue 
turning point (K) approaches the singularity zk ^ ji0 + iP)- 

^'^ This is just the same as the Polya vector field for k,{z,u) represented on the r(z) complex 
plane. The r coordinate is simpler to use in numerical applications 
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Fig. 13: Polya plot of n{r,u)f{r)^^ for /i = 4, d = 8 and u = 1.7e^*T5, 
u = 1.82e^*T6 and u = 2e^^~> , the lines in color are sketched for visual purposes 
and do not represent a faithful integration of the field 




Fig. 14: Same as above but for u = 2e and u = 2e * s 

For \u\ large and some value of arg(u) we have Re{Z{zT, Zk)) = (that is figure 4 
above) and the Stokes lines configuration in the z plane is easiest to picture: 




Also pictured is a schematic representation of the active region. In this case too an 
infinite succession T„ and of images of the turning points is active. 

Given the configuration of turning points we have to apply the connection formula 
at T K Ti Ki T2 K2 ■■ counterclockvk^ise. Proceeding as before v^^e v^^ill have for z in the 
active region between T and K 
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where A is defined as in (|A.2| ). Then by using the connection formula at K we get: 



By repeating the steps for Ti and Ki we obtain: 

so that for z — ^ +00 taking into account the periodic succession of turning points we obtain 

A / Oe^^f^^^ \ 

luu \ smh[iyup/2) J 

with 7 defined by ( |A.6D . 

From this formula it is easy to see that as the two turning points exchange dominance 
the exponential factor in (^i-fe^''^'^^^'^^-*) is just a phase and the location of the quasinormal 
modes can be found by solving 

Z{zt,zk) =i^{l + 2n), n = 0,1,2,- •• 

When u is tilted toward the real axis T and K move in the direction of the red arrows 
and exp{Z{zT, zk)) ^ 1 For u tilted towards the imaginary axis T and K move in the 
direction of the blue arrow and exp{Z{zT, zk)) ^ 1- 

A. 4. kj^O 

For k E R and small there are new solutions to /(r)(l + ^) =0 lying near the 
singularity but the turning points giving the dominant contribution to gujp{z) are the same 
as for k = except for the fact that their position will depend continuously on k. The 
WKB approximations to gujp{z) are the same as before except that the integral defining Z 
has to be evaluated at finite k. 
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■^a ■^min ^max ■^c Z I 

Fig. 15: The left figure is a scliematic plot for the potential ^(z) for k > fee- On 
the right is a simplified scheme of the active region for lra(u) < Re{u) > 



However as described in section 5.2 for finite n when k becomes greater than a critical 
value /cc(/i) the potential on the real z axis changes from being monotonically decreasing 



to the one depicted in fig. |T5| . There is a region between Zmin and Zmax for which the 
potential is an increasing function of z. 

The resulting pattern of Stokes lines for u close to the real axis and between Umin 
and Umax can be quite complex. However if we ignore the complications arising from the 
periodic structure of the turning point in the z plane we can picture the active region for 
< Umin < Re{u) < Umax and u with a small negative imaginary part as in fig. |15. 

This configuration of active turning points gives for the coefficient of e"*'^"^ in the 
WKB expression for g^jpiz) as 2 ^ oo: 

I -j_ Q2iyZ(za,Zc) _|_ ^2iyZ{zt,Zc) = I -\- ^2iyZ{zt,Zc) _|_ ^2iyZ{za,zi,)^ ^j^ 

Notice that with the choice of signs implied by the arrows in fig. H |e2^^(^-^i>)| = 1+ 5 
with 6 while je^'^^C^b'^c) | = A ^ 1. So whenever 

2uZ{za,Zb) ^ in{l + 2n), n = 0,1,2,- •• (A.8) 



( |A.7| ) is zero as long as A6 = 1 An explicit expression for 6 can be given in terms of the 



oscillation period T{za, Zb) between Za and Zb of a classical particle: 



'2uulm(u) 



k{u, k] z) ^dz 
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-i>uIm{u)T{za, Zb) (A.9) 



Therefore we get for the distance of the quasinormal modes from the real u axis: 



Im{u) 



{uuT{za,Zb)) e 



(A.IO) 



As M — > Umax the imaginary part of the quasinormal frequencies increases as A decreases. 
The mode lines diverge from the real axis and for u > Umax they join with those determined 
in the previous subsection as for large u the non-monotonicity of the potential becomes 
less and less important. 

Appendix B. Tortoise coordinate 

In this appendix the analytic structure of the function r(z) will be described in some 
detail and we will present its determination used in appendix A. 
First recall the definition of the tortoise coordinate. 



where /(r) = + 1 - 

The solutions of /(r) = will be denoted in the following by i = 0, 1, .., (i — 1. In 
particular ro is the positive real solution corresponding to the horizon. Also the residues 
at Ti in the integral for z{r) will be denoted as Res{ri) = Then (3o = (3 the inverse 
temperature of the black hole. 

In order to obtain a definite value of z{r) we have to specify a particular contour 
connecting oo to r. Suppose that one particular contour C is chosen and the value zc{r) 
is obtained. All the possible values of z corresponding to the same point will be z{r) = 
zc{r) + 2i Y^iZo where G Z. 

We can give a unique prescription for z{r) by imposing that the contour C cannot 
cross a set of lines starting from the r^. For example we can take these lines as propagating 
from to the radially but this is only one of many possible choices 111. These lines will 
then be branch cuts for the function z(r). 

Some properties of the function z{r) which will be useful in the following are: 

1 The behavior for r ^ oo is z{r) ~ ^ 

2 The behavior for r ~ is z{r) ~ — ^ log(r — r^). 

note that this choice works due to the fact that the sum of the residues is 0. 




(B.l) 
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3 The behavior for r ~ is 2; ~ {d-i)fi ^^^^ where z{0) is any one of its possible 
values. Viceversa r ~ ((d — l)fj,{z — z{0)))'^ . 

From the last item in the list above we see that the function r{z) is not single valued. 
The following is the description of the particular determination of this function that is 
used in appendix A. Notice that for r ~ ro the function r{z) = ro + exp(^^) is periodic 
with period i^. In fact the lines of constant real part for z{r) around tq are topologically 
Si around tq. This is valid only for Re{z) greater than ^ = drP[f(r)~^] where P 
denotes the principal part at the horizon. At this value of Re{z) the line passes at r = 
for z = ^ + i^{n + ^) n G Z. These points are then branch points for the function r{z). 
We will choose the branch cuts propagating from them to have constant imaginary part 
and to go towards Re{z) = —00. By following this prescription for all branch points we 
find a determination of r{z) such that r{z + i^) ~ r{z). With this choice of determination 
for Re{z) > the only branch points are those described: z = ^ + i^{n+ ^) neZ. 
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